arXiv: 1508.02415vl [cond-mat.mtrl-sci] 10 Aug 2015 


Polarization Dependence of Water Adsorption to 
CH3NH3PIJI3 (001) Surfaces 


Nathan Z. Koocher, Diomedes Saldana-Greco, Fenggong Wang, Shi Liu, and 

Andrew M. Rappe* 

The Makineni Theoretical Laboratories, Department of Chemistry, University of Pennsylvania, 

Philadelphia, PA 19104-6323 USA 

E-mail: rappe@sas.upenn.edu 

Abstract 

The instability of organometal halide perovskites when in contact with water is a serious 
challenge to then - feasibility as solar cell materials. Although studies of moisture exposure 
have been conducted, an atomistic understanding of the degradation mechanism is required. 
Toward this goal, we study the interaction of water with the (001) surfaces of CH 3 NH 3 PM 3 
under both low and high water concentrations using density functional theory. We find that 
water adsoiption is heavily influenced by the orientation of the methylammonium cations close 
to the surface. It is demonstrated that, depending on methylammonium orientation, the water 
molecule can infiltrate into the hollow site of the surface and get trapped. Controlling dipole 
orientation via poling or interfacial engineering could thus enhance its moisture stability. We 
do not see a direct reaction between the water and methylammonium molecules. Furthermore, 
calculations with an implicit solvation model indicate that a higher water concentration may 
facilitate degradation. 
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Introduction 


Solar cells based on organometal halide perovskite (OMHP), especially methylammonium lead 
iodide (MAPM3), have had a remarkable increase in efficiency in the past five years.-""- Although 
the photovoltaic mechanism underlying this high power conversion efficiency is not fully under¬ 
stood, the impressive efficiency of MAPM3 is related to its suitable band gap, good carrier transport 
properties, high optical absorption, and long diffusion length.-"— Despite its promising efficiency, 
commercial use of the OMHP-based solar cell is limited in part by its poor stability with respect 
to moisture. 1314 One proposed work-around involves packaging the OMHP with hydrophobic 
materials to enhance stability, 15-16 but this has considerable drawbacks including more compli¬ 
cated device architecture, additional manufacturing requirements, interfacial defects, and necessi¬ 
tating insight into interfacial atomic and electronic structures. I- 7 " 23 Understanding the mechanism 
of degradation is therefore critical to providing materials design principles and engineering strate¬ 
gies for achieving long-term stability. 

The degradation mechanism is currently an open question. Niu and coworkers proposed that 
methylammonium iodide (MAI) and PbU are first formed, with further breakdown of MAI into 
methylamine (CH3NH2) and hydrogen iodide (HI), concluding with the formation of l2(solid) and 
H2(gas) after exposure to oxygen and sunlight. 24-25 Mechanistically, Frost and coworkers sug¬ 
gested that water abstracts a hydrogen from the MA molecule in an acid-base reaction which leads 
to the formation of HI, CH3NH2, and PbU. 26 Further investigation of the initial steps of the degra¬ 
dation process have been performed. Leguy and colleagues suggested that the MAPM3 responds 
differently to moisture depending on the water concentration based on experimental observa¬ 
tions. 27 When exposed to low humidity, MAPM3 forms a transparent monohydrate (MAPbls^O), 
which can be dehydrated back to MAPM3 by raising the temperature. After prolonged exposure 
to water vapor, the monohydrate converts to a dihydrate ((MA)4Pbl6-2H20), which eventually dis¬ 
solves in water, leading to decomposition. 14-27-28 It is also suggested that water could penetrate 
into the perovskite along grain boundaries and that irreversible decomposition occurs when a grain 
boundary has completely converted to the monohydrate. 27 The proposed hydration of grain bound- 
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aries is not fully understood and will benefit from an atomistic examination of how water penetrates 
into the perovskite. Both Christians et al. and Yang et al. found hydrate formation, although the 
exact nature of the hydrate species is ambiguous. 29,30 Both groups observed at least partial recov¬ 
ery of the perovskite when the system was dehydrated. Moreover, Christians et al. found that the 
perovskite decomposes differently in dark and light conditions; 29 the presence of moisture and il¬ 
lumination cause the perovskite to degrade to PbB, but in the dark, Pbl2 is not formed. It should be 
noted that Philippe et al. do not see evidence of a hydrate state in their photoelectron spectroscopy 
data, but attribute it to the instability of the hydrate phase. 31 Although experiments are providing 
an increasingly clear picture of the degradation mechanism, an atomistic-level understanding of 
the process could uncover new methods to stabilize the material. 

The observation of a monohydrate state suggests that the degradation mechanism starts with 
a surface-environment interaction process. Theoretical surface studies of this system so far have 
investigated stable terminations of tetragonal 32 and orthorhombic 33 MAPM3, as well as adsorp¬ 
tion of anisole (a hole-transport material proxy) on ( 001 ) surfaces of pseudocubic MAPM3, 34 and 
the interaction of water with these surfaces with first principles molecular dynamics 35 and density 
functional theory. 36 In this paper, we present an atomistic perspective of water interacting with 
the relevant PM2- and MAI-terminated ( 001 ) surfaces of MAPM3 possessing different polarity. To 
simulate low relative humidity conditions, we introduce one explicit water molecule on the surface 
per primitive surface cell, and to simulate high relative humidity conditions, we use an implicit 
solvation model with and without an explicit water molecule. We show that surfaces with differ¬ 
ent terminations and polarities respond differently to water, leading to different surface bonding 
configurations. 37,38 We also investigate the effect of water penetration into the material. We find 
that the MA dipole orientation strongly affects the surface-water interaction at low coverages and 
postulate that control of the orientation through poling with an electric field or interfacial engineer¬ 
ing between the perovskite material and capping layers could enhance the moisture stability of the 
material. Also, the implicit solvation model results in the elongation of the lattice vector perpen¬ 
dicular to the surface, suggesting that the dielectric response of a higher water concentration aids 
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in material degradation by expanding the lattice. 

The polarization of a material is known to affect the surface adsorption of small molecules. — — 
Due to the permanent dipole moment of the MA cation, there are two orientations that produce 
polarization extrema: all molecules either having their NH3 + -ends pointing toward the surface, 
termed P + , or all molecules have their CH3-ends pointing toward the surface, termed P . The 
top-down and side views of the bare surfaces considered are shown in Fig. 1 . 


Methodology 

The surfaces were modeled using a slab model within the Quantum ESPRESSO DFT package. 43 
The PBE-GGA functional 44 was used as well as norm-conserving, optimized pseudopotentials 
with a plane-wave cutoff of 50 Ry, generated with the OPIUM code. 45-46 The slab model for both 
the MAI- and Pbl2-terminated surfaces had 9 layers, including the water adsorbate, and a vacuum 
of about 15 A between each slab. Geometry optimization calculations involved relaxing the top 
four layers and the adsorbates, until the forces on atoms were less than 0.01 eV/A in all directions, 
while the bottom five layers were fixed to the tetragonal bulk structure. A 4 x 4 x 1 Monkhorst- 
Pack k-point grid was used for the relaxation, while a finer grid was used as needed for the density 
of states calculations. A linear polarization continuum model of water implemented in jDFTx 47,48 
was used to obtain surface structures under liquid water environment. 
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Figure 1 Water adsorption sites on the MAPM 3 (001) surface, (a) Top-down view of the MAI- 
termination, showing exposed MA molecules and associated iodine atoms. Red letters indicate 
studied water adsorption sites, oH + , / 3 ~l + , and y~l + . (b) Side view of the MAI-terminated P + 
surface, where the NH3 + -end is exposed, (c) Side view of the MAI-terminated P surface, where 
the CH3-end of the MA molecule is exposed, (d) Top-down view of the PM2-termination, showing 
coordinated Pb atoms with associated iodine atoms; both are labeled to facilitate discussion. Red 
letters indicate water adsorption sites, cH + , £ 4 , and £ f . (e) Side view of the Pbl2-terminated 
P + surface. The NH3 + -end of the subsurface MA molecule is pointed toward the surface, (f) Side 
view of the Pbl2-terminated P surface. The CH3-end of the subsurface MA molecule is pointed 
toward the surface. The dipole of the MA molecule in each slab is indicated with a red arrow. The 
black lines represent the surface periodicity studied (y^xy^ R 45 °). Grey: Pb, Purple: I, Blue: 
N, Brown: C, White: H. 


Results and Discussion 


We explore different adsorption sites for water on the surfaces by putting the molecule at the sites 
labelled in Fig. 1 and allowing the system to fully relax. As we will discuss later, in some cases, 
different adsorption sites may lead to the same final structure after relaxation. Sites on the MAI- 
terminated surfaces include: on top of an iodine atom, site a; above a hydrogen of the exposed 
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Table 1 Adsorption energies, E a d s in eV, of water on each of the studied sites of MAI- and Pbl 2 - 
terminated surfaces, with the MA dipoles aligned in the P or P arrangement. Note that on the 
PbF-terminated P surface, water molecules starting at <5 and £ ~ relax to C, . Explanation for 
this is given in Figure 2. 


MAI-termination 

Pbl 2 -termination 

P+ 

P~ 

P+ 

P~ 

site 

£ads (eV) 

site 

£ads (eV) 

site 

Eads (eV) 

site 

£ads (eV) 

a + 

-0.49 

a 

-0.18 

5+ 

-0.40 

8~ 

-►c- 

/3 + 

-0.36 

r 

-0.13 

£+ 

-0.38 

£ 


7 + 

-0.48 

7 

-0.28 

c + 

-0.39 

C- 

-0.54 


group of the MA molecule (ammonium, NH 3 + for P + ; methyl, CH 3 for P ~), site /3; and above a 
hydrogen of the lower group of the MA molecule (CH 3 for P + , NH 3 + for P~), site y. The corre¬ 
sponding sites on the P and P + surfaces are denoted with “ and + , respectively and are termed 
a /3~l + , and y + . Sites on the PbF-terminated P surface include: along the Pb-I bond with 
oxygen pointed toward Pb and H pointed toward I, site 8 ; above the hollow site of the surface, 
site e ; and above the Pb, sites . The corresponding sites on the P + surface are <5 + , £ + , and £ + , 
respectively. The adsorption energies of water at each site after relaxation are evaluated as E ads = 
^surface— h 2 o ~ l^hare -surface + F’lboK and are shown in Table 1. Water adsorption is favorable on 
all surfaces, although distinctly different interactions are found depending on both surface termina¬ 
tions and polarities. Adsorption to the MAI-terminated P + surface is more energetically favorable 
than to the corresponding P surface. For PbF-terminated surfaces, however, the water molecule 
prefers to bind to the P surface rather than to the P r surface. This behavior is likely due to the 
competing hydrogen bond interactions between the MA molecules, the Pbl 2 inorganic lattice, and 
the water molecules. For the MAI-terminated surface, water prefers to bind to the P + surface with 
exposed NH 3 " 1 " groups, as a hydrogen bond is formed between water and NH 3 + groups, while no 
hydrogen bond is formed between water and CH 3 groups (P surface). On the Pbl 2 -terminated P 
surface there is only weak interaction between the subsurface methyl groups and the surface PbF 
lattice, so the interaction between water and the surface Pb is stronger than that on the P + surface, 
where the NH 3 1 groups hydrogen bond to the inorganic lattice . 49 Because all starting sites went to 
the same local minimum, we produced a contour plot which shows E a ^ s as a function of position 
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in the (001) plane, seen in Figure 2, indicating that there is only one local minimum around the Pb 
atom. The atomic and electronic structure of the system resulting from various adsorption events 
will now be described. 



a 

Figure 2 Contour plot of water adsorption on the Pbl 2 -terminated P surface of MAPbF,. 

Only one local minimum is seen at the Pb atom, explaining why all sampled structures in Table 1 
resulted in the same adsorption energy. 


Starting with the MAI-terminated surfaces, the lowest-energy configuration of water on the P + 
surface is found from initial site a + shown in Fig. 3a (side view). The water is almost flat in the 
(110) plane with a hydrogen of water (%) interacting with a surface iodine, but is slightly tilted 
to make a hydrogen bond between the oxygen of water (O) and one hydrogen atom (H\) from 
the MA NH 3 + group. To examine the bonding character in more detail, we compute the orbital 
projected density of states (PDOS), seen in Supplementary Information Figure SI. It shows the 
overlap between the O p y orbital and the s orbital of Hn, indicating a hydrogen bond interaction. 
On the MAI-terminated P surface, however, the orientation of the water molecule (as shown in 
Fig. 3b) is very different from that on the P + surface. In this case, relaxation with a water molecule 
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initially at adsorption site a yields a structure where O and one Hw are located directly above 
the surface iodine, with the Hw 2.79 A above the iodine vertically. The PDOS of this structure, 
displayed in Figure SI, shows that the Hw s orbital has weak overlap with the I p orbital and that 
the O p orbitals are split. Comparing the PDOS highlights the weaker bonding of water to the P 
surface. The water-P + surface configuration is 0.31 eV lower in energy than the water-/ J surface 
configuration, probably due to the strong electrostatic repulsion between the oxygen and iodine on 
the P surface without the compensation of the hydrogen bond. 



Figure 3 Water adsorption to the MAI- and Pbl 2 -terminated surfaces with different polari¬ 
ties. Lowest-energy structure of water on the (a) MAI-terminated P + surface, (b) MAI-terminated 
P surface, (c) PbH-terminated P 1 surface, and (d) Pbl 2 -terminated P surface. 


These results suggest that the MAI-terminated P surface is more water resistant than the P + 
surface due to the lack of water reactivity with the CH 3 -end of the MA molecule. Thus, exposing 
the CH 3 to water by applying an electric field 50 or by inducing local, interfacial ordering of the 
dipoles with a capping layer 19 could improve stability. If, however, the water molecule is approx¬ 
imately 2.00 A from the NH 3 + -end of the MA molecule of the P surface (adsorption site y ), 
the MA dipole locally orients to favor the NH 3 + hydrogen bonding with the water, leading to a 
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0.093 eV lower-energy configuration of water on the P surface. Based on previous molecular 
dynamics studies, water travelling from vacuum to this distance is likely spontaneous. 35,36 This 
indicates that the surface is not water resistant in all cases, implying that an electric field might 
have to be applied continuously in the presence of water to improve stability. In general, we see 
no evidence of chemical reactivity between the tbO and the MA molecule, leaving the molecular 
cations intact regardless of orientation or proximity. This is consistent with first principles molec¬ 
ular dynamics results showing that degradation is a structural rather than acid-base effect. 35 We 
cannot rule out an acid-base mechanism, however, because we do not sample the full pH range in 
our calculation. 26,27,29,30 

The lowest-energy water adsorbate configuration on the Pbl 2 -terminated P surface is shown 
in Fig. 3c. This structure is obtained from the relaxation of the water molecule initially on site 
8 (along a Pb-I bond with the oxygen of water pointing toward the Pb). The water molecule 
becomes almost planar in the (001) plane, with the O making a 2.68 A bond with the Pb, while 
one hydrogen from the water interacts with a surface iodine 2.91 A away, as indicated in Fig. 3c. 
On the Pbl 2 -terminated P surface, the water molecule starting on adsorption site 8 1 ends up in 
a similar configuration as that on the P + surface, but is closer to the surface, as shown in Fig. 
3d. The water is oriented almost flat in the (001) plane above the surface, with the Hw atoms 
pointing slightly downward for better bonding with the surface iodine atoms. The O makes a 2.62 
A bond with the surface Pb, while the Hw atoms interact with iodine atoms about 2.80 A away. 
Interestingly, these distances are shorter than in the case of water on the corresponding P + surface 
site, indicating stronger bonding. This message is reinforced by the splitting of the Pb p orbitals in 
the P~ case, with no such splitting in the P + case, as shown in Figure SI. The change in bonding 
character is also reflected in the 0.14 eV larger adsorption energy of water on the P surface than 
on the P + surface (Table 1). 

Experimental observations indicate that water molecules penetrate into the material along grain 
boundaries, converting the structure to the monohydrate state. 27 Therefore, we explore the possi¬ 
bility of water penetrating through the hollow site in the Pbl 2 -terminated surfaces. We model two 
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Figure 4 Energetic and structural description of possible water penetration mechanism on 
the Pbl 2 -terminated surfaces with different polarities, (a) Lowest-energy structure of water 
at the hollow site of the PbL-terminated P + surface, (b) Lowest-energy structure of water under 
the first layer of the PbL-terminated P + surface, (c) Reaction path of water traveling from the 
configuration seen in (a) to that in (b). The forward activation barrier is 0.023 eV, which can be 
overcoming using thermal energy. This emphasizes the facility of water penetration on this surface, 
(d) Lowest-energy structure of water initially placed at the hollow site of the Pbl 2 -terminated P 
surface. The water was repelled from the hollow site into the vacuum above the hollow site during 
relaxation, (e) Lowest-energy structure of water under the first layer of the Pbl 2 -terminated P 
surface, demonstrating that if water can penetrate the top layer, it will tilt the methylammonium 
dipole direction to form a hydrogen bond and become trapped, (f) Reaction path of water traveling 
from the configuration seen in (d) to that in (e). The forward activation barrier is 0.27 eV, which is 
an order of magnitude higher than on the surface of opposite polarity. This emphasizes the water- 
repelling nature of the surface and suggests paths of material stabilization based on poling with an 
electric field or interfacial engineering to create an ordered cation domain. 
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scenarios on both P + and P surfaces: H 2 O in the plane of the first surface layer, and H 2 O inside 
the material (under the first surface layer). Note that this is different from adsorption site £ in Fig. 
1 because £ corresponds to water in vacuum above the hollow site. On the Pbl 2 -terminated P 1 
surface, the water that starts in the plane of the surface stays in-plane, with the O interacting via 
hydrogen bonding with the Hn of the subsurface MA molecule, while the two Hw atoms interact 
with surface iodine ions. This is shown in Figure 4a. The PDOS of this structure displayed in 
Figure S2 shows that the O p x orbital and the Hn 5 -orbital has appreciable overlap. This indicates 
that there is a moderate interaction between the two atoms, weakening octahedral rotations. The 
Pb-I-Pb bond angles are closer to 180°, leading to a larger hollow site. Comparing the top-down 
view of the bare surface (Fig. S3) with that of this structure, it is clear that H 2 O infiltration in¬ 
creases the size of a hollow site, and that the rotation suppression opens up the adjacent hollow 
sites, suggesting that water adsorption on the subsurface has a collective effect on the neighboring 
structure. 

Manually moving the water through the hollow site of the P + surface, such that the water 
gradually approaches the subsurface, yields a trapped, intact molecule (shown in Fig. 4b). The 
PDOS for this structure (Figure S2) shows that the O p orbitals are split, with the p z orbital having 
appreciable overlap with the s orbital of Hn atom in the MA molecule, signaling strong hydrogen 
bonding. To ensure that the structures in Figs. 4a and 4b are distinct minima, we performed a 
Nudged Elastic Band (NEB) calculation with those structures as endpoints in the path, shown in 
Fig. 4c. We found that there is a 0.02 eV activation barrier going from configuration in Fig. 4a to 
that in Fig. 4b, while there is a 0.05 eV barrier for the reverse direction. These barriers suggest that 
water could become trapped beneath the surface layer, with a facile equilibrium at room T. 

The water was also placed in the P surface plane at the hollow site. As expected based on the 
previously observed hydrophobicity of the methyl group, instead of staying in-plane, the water was 
repelled from the surface. As seen in Fig. 4d, the optimized water orientation has the O pointed 
away from the surface and the two Hw atoms interacting with the surface iodine (Figure S2). 
When the water is between the first and second layers in the Pbl 2 -terminated P system, the water 
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molecule is trapped in part due to hydrogen bond formation (Figure S2) with the NH 3 + -end of 
the methylammonium, which is achieved through a change in orientation of the methylammonium 
cation. We also calculated the activation energy barrier using NEB for this pair of structures (Figs. 
4d and 4e) and found that the forward barrier is 0.27 eV and the reverse barrier is 0.57 eV, seen in 
Fig. 4f. As expected, the forward barrier for this surface is higher than that for the corresponding 
P + system and once the molecule penetrates the hollow site, it cannot leave. Interestingly, the 
forward barrier is more than ten times the thermal energy, further indicating that the methyl end of 
the methylammonium can repel water. 

Contrasting these surfaces, we see that the oxygen binds more strongly with the lead on the P 
surface than on the P + surface and interacts more favorably with the ammonium end of the MA 
molecule than with the methyl end regardless of the surface termination. On the P surface, the 
water molecule is repelled from the hollow site, while on the P + surface, the water molecule is 
stable in-plane with the top atomic layer, probably due to the strong hydrogen bond between the 
water and the NH 3 + of the MA molecule. Furthermore, both surfaces can trap water molecules in 
the subsurface, and such trapping disrupts the Pb-I-Pb angles of the surface layer. These results 
suggest a possible scheme of how water interacts with the surface and possible methods to stabilize 
the material. On approach to the surface, the water would first get physisorbed in a shallow energy 
well, forming a hydrate. Then, environmental factors such as temperature and humidity could 
enable the FFO to move to the surface, where it could get trapped at the interfacial plane of the 
hollow site on the P + surface ,- 5 but still be repelled from the hollow site by the P surface. 
If the water molecule penetrates the hollow site of the surface (of both polaritities), water and 
the Pbl 2 inorganic lattice compete to form a hydrogen bond with MA, weakening the stability 
of MAPM 3 . Most importantly, this suggests that the stability of MAPbl 3 -solar cells could be 
enhanced if they were poled using an electric field or locally ordered via an interfacial substrate 
such that the methyl ends of the MA molecules were pointed toward the surface. While our ki netic 
barrier results demonstrate that water incursion into the material is slowed, without knowing other 
steps leading to chemical reaction, it is difficult to assess the actual extent of the protection. More 
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work is necessary to reveal the effect of electric field poling and interfacial engineering on material 
stability. 


Table 2 Reaction energetics (eV) of dissociated and oxidative water on Pbl 2 -terminated surfaces 
with different polarities 




P+ 

P~ 

Pbl 2 surl 

Pbl 2 surf 

+ H 2 0 -► HO-PbI 2 surt + {Hi 
+ H 2 0 0=PbI 2 surf + H 2 

1.46 

3.02 

1.52 

4.35 


In the interest of exploring beyond molecular adsorption of water, we constructed dissociated 
and oxidative water configurations on the PbR-terminated P+ and P surfaces. The optimized 
structures can be seen in Figure S4. The first case simulates PRO dissociation by placing OH and 
H+ apart from each other. The optimized structure for both P 1 and P~, however, showed water 
convert from an initial dissociated state to a molecular form exhibiting similar bonding features to 
the previously described molecular adsorption configurations. Similarly to the structures described 
before, the oxygen was closer to Pb on the P than on the P+ surface. We also studied a lone hy¬ 
droxyl bonded to Pb to model oxidative adsorption. The hydroxyl causes the Pb to rise up out of 
the surface. Finally, we investigated a single oxygen atom for stronger oxidative adsorption. On 
the P + surface, the oxygen atom infiltrated the PbF layer, bonding between the Pb and I. On the 
P surface, the oxygen was incorporated into the surface plane, disrupting the periodicity of the 
PbH-tcrminatcd surface indicating the formation of PbF units. These oxidative adsorption con¬ 
figurations were higher in energy than molecular adsorption configurations (see Table 2), making 
them unlikely under low humidity conditions. 

Table 3 Elongation of top interlayer spacing, / SU rf-subsurf, of Pbl 2 - and MAI-terminated surfaces 
and shortening of the vertical distance between FFO and PbR-terminated surface, /?n 2 o-Pb- with 
different polarities computed without and with polarizable continuum model (PCM) as an implicit 
solvation model for water. 



Without PCM 

PCM 

Change after PCM is applied 

Pbl 2 -terminated P+, / SU rf-subsurf (A) 

6.35 

6.48 

2.0% 

MAI-terminated P+, /surf— subsurf (A) 

6.83 

7.08 

3.7 % 

PbH-terminated P+, £>HiO-Pb (A) 

2.68 

2.58 

-3.7 % 

Pbl 2 -terminated P~, /^HiO-Pb (A) 

2.57 

2.50 

-2.7 % 
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To approximate the effect of high humidity conditions on the MAPM 3 surface, which is compu¬ 
tationally expensive to treat with explicit water, we employ a polarizable continuum model (PCM) 
as the implicit solvation model on the MAI- and Pbl 2 -terminated bare P + and P surfaces, as well 
as on the Pbl 2 -terminated P + and P surfaces with one water molecule added explicitly. These 
calculations were performed with the Joint Density Functional Theory package.—— The bare sur¬ 
faces show an elongation in the [ 001 ] (out-of-plane) lattice constant upon inclusion of solvation, 
as shown in Table 3. The Pbl 2 -terminated P 1 and P surfaces with explicitly adsorbed water show 
a shorter O-Pb bond length, indicating a stronger bond. The expansion normal to the surface when 
solvation is included suggests that higher water concentration weakens the interlayer bonding of 
the OMHP 

Conclusion 

In this work, we present an atomistic view of water interacting with MAPM 3 (001) surfaces. We 
find that water favorably adsorbs on all the studied sites in both MAI- and PM 2 -terminated surfaces, 
supporting the existence of the experimentally observed hydrate state as a potential initial step 
of degradation of the material. Our results do not show a direct reaction between the FPO and 
MA molecule, which has been proposed to be an initial step of degradation . 26 Although this is 
consistent with previous work , 35 we cannot completely rule out an acid-base reaction. The specific 
bonding characteristics of water with the OMHP are termination and polarity dependent, with the 
MAI-terminated P + surface having more favorable water adsorption than the corresponding P 
surface, while the PbH-terminated P surface binds water more strongly than the corresponding 
P + surface. These interactions arise from the hydrogen bond of the ammonium group of the MA 
molecule with the inorganic lattice and the water molecules, while the hydrophobic methyl group 
does not react with the water. Harnessing the water-repelling nature of the methyl group of MA by 
poling the surface during device operation or creating a local domain through interface engineering 
could thus enhance the stability against degradation, although this stabilization strategy requires 
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further investigation. 

We demonstrate that physically introducing a water molecule below the top atomic layer of 
Pbl 2 -terminated surfaces pushes up the surface layer and affects the adjacent hollow sites, provid¬ 
ing important insight into how the material degradation could propagate. The presence of more 
oxidizing adsorbates, though unlikely to form in low humidity conditions, has a greater effect 
on the surface structure. Finally, calculations with an implicit solvation model of water provide 
a more general picture of water effects and show that the lattice elongates perpendicular to the 
surface in the presence of FPO. This implies that a higher water concentration may weaken the 
interlayer bonding strength of the hybrid perovskite, facilitating the dissociation of MAPM 3 into 
its constituents. 
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